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ABSTRACT 

The fragmentation process in collapsing clouds with various metallicities is studied 
using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by 
three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For 
different combinations of these parameters, we calculate 480 models in total and study 
cloud evolution, fragmentation conditions, orbital separation and binary frequency. For 
the cloud to fragment during collapse, the initial angular momentum must be higher 
than a threshold value, which decreases with decreasing metallicity. Although the exact 
fragmentation conditions depend also on the initial cloud shape, this dependence is 
only modest. Our results indicate a higher binary frequency in lower-metallicity gas. 
In particular, with the same median rotation parameter as in the solar neighbourhood, 
a majority of stars are born as members of binary/multiple systems for < 10~ Zq. 
With initial mass < 0.1 M©, if fragments are ejected in embryo from the host clouds 
by multi-body interaction, they evolve to substcllar-mass objects. This provides a 
formation channel for low-mass stars in zero- or low- metallicity environments. 
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1 INTRODUCTION 

Observations in the sola r neighbourhood have shown a high 
binary frequency (lAbtl 1 19831 : iDuquennov fc Mayor] Il99ll: 
Fischer fc Marcvlll992l ; see also the review of lGoodwin et all 
20071 ). Because the binary frequency in star-forming regions 
is even higher than that in the field (|Mathiedll994 ). a ma- 
jority of stars is believed to be born as binaries. However, 
the observable star-forming regions are limited in the solar 
neighbourhood and thus at metallicity comparable to the 
solar abundance (i.e., Z = Zq). We do not know the binary 
frequency in lower (Z < Zq) or zero (Z — 0) metallicity 
environments. 

Binary stars have important roles in many astrophysical 
contexts. They are considered to be related to a pollution of 
extremely metal-poor stars (Suda et al. 2004; Lucatello et al. 
2005; Komiya et al. 2006), progenitors of gamma-ray bursts 
(Bromm & Loeb 2006) and Type la supernovae (Nomoto, 
Thielemann, & Yokoi 1984), the evolution of globular clus- 



* E-mail: machidam@scphys.kyoto-u.ac.jp (MNM); 

omukaiOth. nao.ac.jp (KO); matsuOi. hosei. ac.jp (TM); 
inutsuka@tap.scphys.kyoto-u.ac.jp(SI) 



ters (Sugimoto & Bettwieser 1983; Makino 1996) and the 
source of gravitational waves (e.g., Seto 2002). 

The formation process of binary (or multiple) systems 
in collapsing clouds with Z = Zq has been studied by 
m eans of three-dimensional hyd rodynamics (see the reviews 
of iBodenheimer et alllioool and iGoodwin et al.ll2007l ) . They 
showed that fragmentation occurs after the cloud becomes 
optically thick at n > 10 11 cm~ 3 , leading to formation of 
binary or multiple stellar systems. Whether fragmentation 
occurs depends sensitively on the initial co ndition, e.g. the 
initia l angular velocity and cloud shapes (jGoodwin et all 
l2007t ). Thus, to determine the fragmentation condition and 
binary frequency, we need to calculate cloud evolution for a 
number of models. 

On the other hand, star formation in primordial 
clouds has also been inv estigated by cosmological three- 
dimensional simulations (|Abel et all 120021; iBromm et all 
120021 : lYoshida et a"ill2006l) . Recently, lYoshida et all (|2008l ) 
succeeded in calculating cloud evolution up to protostar for- 
mation (n ~ 10 21 cm" 3 ) and demonstrated that fragmen- 
tation does not occur. However, cosmological simulations 
tend to calculate only a single clump with the lowest angu- 
lar momentum, which collapses first. This setting prevents 
the cloud from fragmenting, because binaries form in clouds 
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with a high angular momentum. Therefore, if we consider 
not only the clouds collapsing first but also those forming 
at a later time, bina ries may form within them. In fact, 
ISaigo et all (12004 ) and lMachida et all (|2008aT ) calculated the 
collapse of primordial clouds with different degrees of initial 
rotation, and sh owed the pos s ibility of fragmentation and bi- 
nary formation. IClark et al.l (|2008h demonstrated fragmen- 
tation of zero- and low-metallicity turbulent clouds by three- 
dimensional hydrodynamics and showed that the fragments 
acquire angular momentum from turbulence. 

The cloud evolution and fragmentation condition only 
for pri m ordial clo u ds wer e investigated in iMachida et al.l 
(|2008tJ ). iMachidal (|2008br > investigated the evolution of 
clouds with various metallicities and showed that the binary 
frequency increases as the metallicity decreases. However, in 
this study, almost spherica l cloud s are adopted as the initial 
states. As iGoodwin et all (|2007T l pointed out, slight differ- 
ences in the initial conditions may affect cloud evolution and 
its fragmentation process. In the present paper, we investi- 
gate cloud evolution with var ious metallicities for a wider 
range of initial conditions than lMachidal l|2008bl ) . Especially, 
to investigate the effect of the cloud shape on fragmentation, 
in some models, we adopt a strongly distorted structure as 
the initial state. In total, we have calculated 480 models with 
combinations of three parameters: the cloud metallicity (Z), 
initial rotation energy (/3o), and initial amplitude of the non- 
axisymmetric perturbation (A^,) corresponding to the initial 
cloud shape. The large number of models makes it possible 
to discuss the fragmentation condition and epoch, separation 
between fragments for clouds with different metallicities. 

This paper is organized as follows. The basics of our 
model are summarized in §2. The thermal evolution of clouds 
with different metallicities is shown in §3. Then, we describe 
the dependence on initial cloud parameters in §4. In §5, we 
discuss some uncertainties and possible fates of the frag- 
ments. Finally, we summarize our results in §6. 



has 7K. As critical Bonnor-Ebert spheres are assumed as 
the initial state, the radii of the initial spheres are differ- 
ent among models with different metallicities (or different 
initial temperatures): the radius for models with Z = 
is 5.5 x 10 5 AU, while that for models with Z = Zq is 
1.2 x 10 AU. We confirmed that initial differences in radius 
do not greatly affect subsequent cloud evolution, in which 
different cloud sizes but the same metallicity are adopted as 
the initial state. 

To induce fragmentation, the non-axisymmetric density 
perturbation of the m = 2 mode (i.e. the bar mode) 8p is 
added to the initial cloud: 

_ f J 4^(r- c /r- p ) 2 cos20 for r c < r p , . . 

1 A^, cos 20 for r c > r p , 

where A^, and r p are the amplitude and effective radius of 
the density perturbation, respectively, and r c is the cylin- 
drical radius on the equatorial plane. We adopt A$ as 
A^ = 10 -4 - 0.4, and r p as 0.1 R c and R c . Models with 
Tp = 0.1 R c have a larger non-axisymmetric perturbation 
than those with r p = R c . The initial cloud has a strongly dis- 
torted shape in models with large A^, and R c = 0.1. For con- 
venience, we describe the amplitude of the non-axisymmetric 
perturbation as A$, in which a capital 'L' is added after the 
value of the amplitude only when r p = 0.1-R C is adopted. For 
example, the model with A^ = 0.1L has parameters of (A^, 
Tp) = (0.1, 0.1i?c), while that with A^ = 0.1 has parameters 
of (Ap, r p ) = (0.1, Rc). 

The models are characterized by three parameters: the 
cloud metallicity Z, initial rotation energy /3o and degree 
of non-axisymmetric perturbation A^,. The initial rotation 
energy means the ratio of rotation to gravitational energy 
in the initial core. The values used for these parameters are 
Z = 0-lZ Q , A, = 10 -6 - 10~\ and A^ = 10 -3 - 0.3 and 
0.1L— 0.4L. Combining three parameters (Z, /3q, and A&), 
we calculated 480 models. 



2 MODEL SETTINGS 

We solve the equations of hydrodynamics adopting a nested 
grid method to ensure the Jeans condition of the collapsing 
cloud. Each level of a rectangular grid has the same number 
of cells of 128 x 128 x 8 (for details, see Machida et al. 2005a; 
2006a). The energy equation is not solved: instead, we adopt 
the barotropic relations for pressure as a function of density 
calculated with one-zone models by Omukai et al. (2005). 
Figure [1] shows these relations for different metallicities as a 
function of the number density. 

As an initial condition, we take a spherical cloud with 
density 1.4 times higher than hydrostatic equilibrium with 
external pr essure (i.e. the so-ca l led c ritical Bonnor-Ebert 
sphere; see iBonnor! Il956l ; lEbertlll955l ). The initial central 
density is taken as n c ,o = 1.4 x 10 4 cm" 3 . The critical radius 
of the Bonnor-Ebert sphere is R c = 6.45 c s /[47tG'psb(0)] 1 / 2 . 
Outside this radius, a uniform gas density of ube(Rc) = 
10 3 cm -3 is assumed. Each cloud rotates rigidly (fio) around 
the 2-axis. The initial temperatures, which are derived from 
the one-zone model, are different in clouds with differ- 
ent metallicities (see Fig. [TJ. For example, a cloud with 
Z = (primordial composition) has an initial temperature 
of 230 K, while a cloud with Z = Zq (solar composition) 



3 OVERVIEW OF THERMAL EVOLUTION 
AND ITS CONSEQUENCES FOR 
PRESTELLAR COLLAPSE 

Before discussing fragmentation process, we briefly summa- 
rize the outline of thermal evolution with different metallic- 
ities (Fig. [TJ , and see its consequence on dynamical evolu- 
tion of clouds with very small rotation energy /3o = 10~ 6 . 
For those models, the radial density and velocity profiles are 
shown in Figures [2] and [3] at some different epochs from the 
initial state (~ 10 4 cm -3 ) to after the protostar formation 
(> 10 21 cm- 3 ). 

The thermal evolution of the clouds with Z — 10~ 6 Z@ 
is almost the same as the Z = case. For the clouds with 
Z = and Z = 10 -6 Zq, the density and velocity distribu- 
tions are almost the same. The temperature remains rather 
high (~ 200 K at ~ 10 4 cm -3 ) and increases gradually with 
the effective ratio of specific heat 7 = dP/dp w 1.1. The 
density gradient in the outer envelope is well approximated 
by oc r~ 2,2 (Fig. [2}. Even though the cloud becomes optically 
thick to the H2 collision-induced absorption at ~ 10 16 cm~ 3 , 
effective cooling by the H2 dissociation prevents the temper- 
ature from rising adiabatically until ~ 10 20 cm (Fig. [1]), 
where the dissociation is almost completed. Subsequently, 
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the growing pressure halts the gravitational collapse at the 
center and the protostar forms, which is surrounded by ac- 
cretion shock (nc £ 10 20 cm" 3 at r ~ 10 12 cm" 3 ; Fig. EJ) 
with a radius and mass of ~ 1 Rq and 10 _3 Mq, respec- 
tively. 

With metal enrichment Z ^ 10~°Zq, a striking dif- 
ference emerges in thermal and dynamical evolution. The 
clouds becomes optically thick to the dust in a density 
around 10 11 - 10 15 cm , which depends on metallicity in a 
way that those with lower Z becomes adiabatic at a higher 
density, and experience a transient adiabatic phase up to 
the onset of H2 dissociation at ~ 2000K (first adiabatic 
phase, hereafter). Thus, during the first adiabatic phase, 
a transient hydrostatic core, called the first core, appears 
(indicated by the dotted arrow in Fig. [2}. The first core 
has a size of ~ 1 AU, and mass of ~ 0.01 Mq at its for- 
mation. The central temperature of the first core increases 
with the mass of the first core, which grows by accretion. 
At ~ 2000K (~ 10 16 cm" 3 ), the H 2 begins to dissociate and 
another episode of dynamical collapse begins. This is called 
the second collapse phase, in contrast with the first collapse 
phase before the first adiabatic phase. The evolution there- 
after coincides with the zero-metallicity case. Following the 
end of dissociation (~ 10 20 cm -3 ), the temperature increases 
adiabatically again (the second adiabatic phase) and a sec- 
ond hydrostatic-core (the protostar) forms with a similar 
size as that in zero-metallicity case. Note that two hydro- 
static cores (the first core and the protostar) appear in a 
nested manner, i.e., the protostar forms inside the first core. 
Correspondingly, two shocks are seen for Z ^ 10 _5 Zq in 
Figures [2] and [3] This existence of two adiabatic phases, one 
at n c ~ 10 11 cm -3 and the other 10 20 cm -3 , during the col- 
lapse has an important consequence in fragmentation nature 
of metal-enriched clouds. 

In the model with Z — 10~ 4 Zq, three shocks (i.e., 
three velocity peaks) appear with an extra weak shock 
at n c — 10 9 cm" 3 by the rapid H2-formation heating at 
n c ~ 10 8 cm -3 (see, Fig. [TJ. However, because of the very 
short duration, no clear adiabatic core is observed at this 
epoch. 



4 DEPENDENCE ON INITIAL CLOUD 
PARAMETERS 

In this section, we see the evolution of clouds with differ- 
ent values of rotation energy /3o and non-axisymmetric per- 
turbation A^, (i.e., cloud shape) at the initial states. For a 
given metallicity, the number of combination (f3o,A$) we 
calculated amounts to 80. With eight different metallicities 
Z = — Zq , 480 models are calculated in total. The density 
distributions on the equatorial plane at the end of calcu- 
lation are presented for different combinations of (f3o,A^), 
which are indicated by the ordinate and abscissa, in Fig- 
ures 2]-[H for a given metallicity. The density perturbation 
is imposed at the scale of the critical radius of the BE sphere 
R c for models shown in the first-seventh columns, while at 
0.1-Rc for those in the eighth-tenth columns (indicated by 
"L"; see, The latter corresponds to models of larger 
non-axisymmetric perturbation (i.e. more distorted). 



4.1 Classification of the Fate of Clouds 

The white-and-black-dotted lines in Figures |4]-[TT] mark the 
border between models in which the clouds fragment or not. 
Each class of models can be classified into two sub-classes, 
which are indicated by the background colours: the fragmen- 
tation models include "fragmentation" (red) and "merger" 
(violet) models, while non-fragmentation ones include "non- 
fragmentation" (blue) and "stable-core" (grey) models. 

Fragments survive without merger in "fragmentation" 
models, while they merge to form a single core before the 
end of calculation in "merger" models. The protostars form 
without fragmentation in the "non-fragmentation" models, 
while, in the "stable-core" models, the first core remains sta- 
ble for 3> Wts, where tg is the local free-fall timescale at 
the centre. Due to our CPU-time limitation, we failed to fol- 
low further cloud evolution for most "stable-core" models. 
We expect that the first core collapses eventually to form a 
protostar after angular momentum transfer in a long-term 
calculation. In fact, for some "stable-core" models, however, 
we succeeded in following the second collapse and protostar 
formation and confirmed no fragmentation. Thus, in this pa- 
per, we regard the "stable-core" models among the cases of 
non-fragmentation when we discuss the fragmentation con- 
dition. 



4.2 Fragmentation Frequency for Different 
Metallicities 

Investigation of Figures 111 - 1 1 1 1 indicates : 

(i) the first parameter to determine whether the a ro- 
tating spherical cloud fragments or not is the initial rota- 
tion parameter /3o- Although the behavior of the boundary 
between the two cases is very complicated, the threshold 
value of /3o for fragmentation increases with metallicity, be- 
ing p = 10" 5 for Z = while (3 a = 10" 3 for Z = Z Q . As a 
result, fragmentation is observed in more models for lower 
metallicity; 23/60 models for Z = Zq, while 47/60 models 
for Z = 0. 

(ii) clouds with large non-axisymmetric perturbations 
develop spiral patterns, which transfer angular momenta by 
gravitational torque and effectively reduce the rotation pa- 
rameters j3o- For example, upward shift of the fragmentation 
boundary towards larger A$ is seen in the Z = 10" 3 - lZ e 
cases. Similarly, the spirals remove angular momenta from 
the rotation-supported first cores, which appears for 0.1, and 
IZq (Figs. [4] and [5}, enabling further collapse to the proto- 
star in large A^, cases. Since the non-axisymmetric pertur- 
bation can grow until the higher density phase for a cloud 
with lower metallicity, fragmentation tends to occur through 
the bar configuration. 



4.3 Spin-up and Fragmentation of Clouds during 
the Collapse 

As seen in Figures [3] - 1111 the fragmentation condition de- 
pends strongly on the initial rotation rate (/3o) but only 
weakly on the initial cloud shape (A$). So far, we used the 
rotation energy of the whole cloud for easier comparison with 
observations. It is, however, known that a more important 
index for cloud evolution and fragmentation is the angular 
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velocity normalized by the freefall timescale (hereafter the 
normalized angular velocity), 

u c = ^p=, (2) 
V47T G p c 

where Q c and p c axe the angular velocity and density 
at the centre of th e cloud (Matsumoto fc Hanawa 2003, 
iMachida et~ai]|2004 iMachida et~aH l2005bl l . For a uniform 
sphere with rigid rotation, the normalized angular velocity 
is related to the rotation energy as 



/3c = 



Q. 2 C R 3 



(3) 



3GM 4nGp c 

where R and M are the radius and mass of the sphere, re- 
spectively. Along with the collapse, the clouds spin up and 
uj c increases without significant angular momentum trans- 
fer. In the case of the primordial (Z = 0) clouds, Machida 
et al. (2008b) found that if the normalized angular momen- 
tum reaches the critical value 0.2 - 0.3 (the gray zone in 
Fig. I12[) and thus a thin disk forms owing to centrifugal 
force before the protostar formation, the disk subsequently 
fragments into binary or multiples. On the other hand, if 
the cloud becomes adiabatic and protostar forms before this 
condition is met, the cloud recovers a spherical shape and 
does not fragment thereafter. 

Here we show that, for higher metallicity, appearance of 
the first core at lower density put more stringent condition 
on lu c for fragmentation. Figure [T^] shows evolution of the 
normalized angular velocity u>c against the central number 
density for all cases. Initially, while the centrifugal force is 
not important, the clouds collapse almost spherically and ui 



(also, f3 c ) increases in proportion to oc 



,V6 



la 1/3 

(Pc oc ric , re- 



spectively). In clouds with rapid initial rotation, u> reaches 
the critical value uif Yag = 0.2 — 0.3 and the clouds frag- 
ment after some more contraction. For Z ^ W~ 6 Zq, the 
critical value must be attained somewhat before the proto- 
star formation (~ 10 17 cm -3 ) for Z ^ IQ~ 6 Z@, as found by 
Machida et al. (2008b), while for Z > 1O~ 5 Z this must 
be before the first core formation (indicated by the vertical 
lines in Fig. I12p . Since the first cores form at lower density 
than the protostars, the clouds have a shorter density range 
to amplify uj in the cases of Z ^ 10~ 5 Zq than the primor- 
dial case. In addition, first core appears earlier for higher 
metallicity. Therefore, an initial higher rotation is required 
to cause fragmentation for clouds with higher metallicity. 

Using the relation uj c oc for spherical collapse, the 
condition for the normalized angular momentum is amplified 
from the initial value ujq at density riini to the critical value 
Wfrag before the density n a di where the gas becomes adiabatic 
is 

UlO > ^frag • (4) 

V n ini ) 

This can be translated to the condition on initial rotation 
parameter Pq: 

/30>Ara g ( — (5) 

where /3 frag = cj f 2 ag = 0.04 - 0.09. 

The adiabatic density n a di and critical rotation param- 
eter /3o,crit given by the right-hand side of equation (JSJ) 
are listed in Table [T] for different values of Z. Here, we 
used a conservative value of /3f ra g = 0.1 is adopted. For 



Z ^ 1Q~ 5 Zq, n a di is the density where the clouds become 
optically thick to dust grains and the first cores form. For 
Z ^ W~ 6 Zq, where no first adiabatic phase is present, 
after reaching /fcag, the clouds need some more density 
interval before the protostar for fragmentation and thus 
i a di = 10 17 cm -3 is adopted as the critical density (see 
Machida et al. 2008a). These values for /3o,crit are in con- 
cordance with the results presented in Figures l4l- II 1 1 



4.4 Fragmentation Epochs and Separations 

Figure [13] shows the separations at the fragmentation epoch 
for all fragmentation models. In each panel, the solid line 
indicates the Jeans length, while the dashed vertical line 
indicates the beginning of the first adiabatic phase. Note 
that the separations are about 10 - 100 times the Jeans 
length since the radial size of the cloud is 10 - 100 times the 
Jeans scale at fragmentation owing to rotation. With the 
vertical scale comparable to the Jeans length, the height-to- 
radius ratio is about I/I0 at this epoch. 

For Z = Zq, fragmentation events cluster near the up- 
per left corner n c = I0 11 — 10 14 cm" 3 , i.e., just after first 
core formation. For lower metallicity, this distribution ex- 
tends toward the lower right because the first core forms 
later and thus fragmentation occurs at a higher density with 
narrower separation for lower metallicity. For Z = 1Q~ Z@ 
as well as W~ 2 Zq, one fragmentation event is observed at 
very high density n c ~ 1 20 cm -3 while the rest locate in 
10 10 cm -3 < n c < 10 14 cm -3 , which are due to the first-core 
formation. In some models of 10~ 4 Zq, fragmentation is ob- 
served also before the first-core formation. The abrupt tem- 
perature rise at n c — 1 8 cm -3 in the 1Q~ 4 Z@ clouds (Figs. [5] 
and [3]) temporarily slows the collapse, enabling fragmenta- 
tion at such an early phase. Without the first-core forma- 
tion, no clear fragmentation epoch exists for Z ^ 10~ 6 Zq. 
Fragmentation takes place whenever a disk-like configura- 
tion appears due to rotation and distributes in a wide range 
of 1 10 cm" 3 < n c < 10 21 cm" 3 (Machida et al. 2008a) al- 
though some events at n c > 10 20 cm~ 3 might be caused by 

5 Z G 



the protostar formation, as in the models with Z ^ 10 
at first-core formation. 

Figure [T3] shows the frequency distribution of further- 
most separations between fragments at fragmentation This 
shows that binaries with higher metallicity have a wider 
separation. Those with Z = Zq have separations 10AU< 
r scp < 1000 AU, while at Z = 0, the range is r scp < 10 AU. 
The separations are distributed in a wide range of 0.1 AU < 



r scp < 1000 AU for models with Z = 10" 



lO'^Z, 



0- 



4.5 Non-axisymmetric Perturbation and 
Fragmentation Mode 

In a few cases in Figures [3] - 111! the clouds fail to frag- 
ment even if the condition (eq. [4j is fulfilled. With moderate 
non-axisymmetry, the cloud can avoid fragmentation even if 
condition (eq. Q is fulfilled owing to angular momentum is 
effectively transferred by the spiral arms. Interestingly, with 
even larger A$ , the clouds tend to fragment through bar con- 
figurations instead of ring configurations realized for frag- 
mentation in small A$ cases. Lower metallicity clouds tend 
to fragment at higher density. This allows non-axisymmetric 
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perturbations to grow for a longer time. Thus, fragmenta- 
tion through the bar is more common for lower metallicity. 
For example, for Z = 0, 33 out of total 47 fragmentations 
are through the bar mode. 



5 DISCUSSION 

As shown in H4.3I the fragmentation condition is that the 
initial rotation energy /3o must be larger than the threshold 
value, which increases with the metallicity. Thus, with lower 
metallicities, clouds with smaller rotation energy can frag- 
ment. Observations in the solar neighbourhood have shown 
that molecular clouds (with solar metallicity) have rotation 
energies in the range of 10~ 4 < ftp < 0.07 with an aver - 
age value of /3 = 0.02 (|Goodman et al.lll993l ; ICasellill2002f ). 
On the other hand, we cannot observe the initial rotation 
energies of clouds with lower or zero metallicity. Cosmolog- 
ical simulations have shown that host clou ds of star for- 
mation have rotation en ergies of (3 ~ 0.1 (|Bromm et alj 
l2002l ; lYoshida et alj2006h . which is slightly larger than those 
in molecular clouds. Thus, we expect a higher binary fre- 
quency in a lower-metallicity environment (or in the early 
universe), roughly assuming that host clouds have almost 
the same distribution of angular momentum (or rotation 
energy) as observations in the solar neighbourhood. If the 
lower limit of the rotation energy is /3o — 10~ 4 , as in the 
solar neighbourhood, fragmentation always occurs in clouds 
with Z ^ 10~ 4 Zq, because their fragmentation condition 
is /3o > 5 x 10 -4 (see, Table 1). As a result, the binary fre- 
quency increases as the metallicity decreases, and most stars 
are born as binary or multiple systems in low-metallicity en- 
vironments. 

In our calculation, fragments have masses of 0.1 — 
10~ 3 Mq at fragmentation. This mass range corresponds to 
that of brown dwarfs or low-mass stars. The mass of frag- 
ments increases in the subsequent accretion phase. Since the 
gas supply (or gas accretion) stops at ejection, ejected frag- 
ments with mass of 10~ 3 — 0.1 Mq are expected to evolve 
to form a brown dwarf or low-mass star. Thus, our result 
indicates that even in extremely low- or zero-metallicity en- 
vironments, substellar-mass objects can form. On the other 
hand, fragments remaining in the host cloud are expected to 
increase their mass by accretion and grow to massive stars 
(or massive star binaries). As a results, we expect that a sin- 
gle massive star (or a massive binary system) and multiple 
low-mass stars can be formed in a single host cloud. When 
multiple fragments appear like a star cluster, as shown in the 
upper left panels of Figures HI - 1111 low-metallicity stars are 
expected to have two peaks at 10 -3 — 0.1 Mq and 100 Mq in 
mass distribution. Note that if fragments stay in the cloud 
for a while and are ejected after acquiring some mass, a star 
with mass M ^ 0.1 Mq can form. 

Subsolar-mass stars can form even in low- (or zero-) 
metallicity environments by the ejection of fragments. Re- 
cently, metal-poor stars with Z < 10~ 5 Zq are observe d 
(|Christlieb et alj|2002l : iFrebel et al.ll2005l ; lAoki et alj|2006t ). 
Forming in the early universe, they should have a mass of 
<J 0.9 Mq. Such low-mass metal-poor stars can form by 
the fragmentation and ejection mechanism, as described 
above. In addition, our result allows even the existence of 
zero-metallicity stars with sub-stellar mass. Our calculation 



showed that the binary frequency for zero-metallicity stars 
is as high as for those with Z = 10~ 5 and W~ 6 Zq. Thus, we 
expect observation of zero-metallicity stars in the future. In 
addition, as in the cases of (A$, /3) — (0.001, 0.1) in Figure[6] 
the binary system can be ej ected from t he ho st cloud when 
multiple fragments appear. ISuda et al.l (|2004l ) claimed that 
stars with Z < 10~ 4 Zq are in binaries, because their surface 
abundance ratios are well explained by mass transfer from 
an (unobserved) AGB companion. The ejection mechanism 
may explain the origin of such low- or even zero-metallicity 
binaries. 

In this study, we followed the evolution of fragments and 
protostars only in their early phase. To determine the final 
mass and fate of fragments, calculations in the subsequent 
accretion phase are required. 



6 SUMMARY 

In this research, we studied the binary /multiple formation 
process during the collapse of rotating clouds with various 
metallicities. We found the following results: 

(i) A transient hydrostatic core ( "first core" ) appears dur- 
ing the prestellar collapse for metallicity ^ 10 -5 Zq at higher 
density and with smaller size toward lower metallicity. For 
Z ^ 10~ 6 Zq, a protostar directly forms without any tran- 
sient core. 

(ii) When rotation reaches a fraction of Keplerian value 
before formation of the adiabatic core (i.e., the first core for 
^ W~ 5 Zq or the protostar for ^ 10~ 6 Zq), the cloud de- 
velops disk-like structure and subsequently fragments into 
binary/multiples. This condition can be translated to the 
initial rotation parameter /3o exceeding a threshold value 
A),crit- Since lower-metallicity clouds have longer density in- 
terval to spin up, i.e. higher density for the adiabatic-core 
formation, they have smaller threshold value. For example, 
A),crit ~ 10~ 3 (10~ 5 ) for Z = Zq (Z = 0, respectively). 

(iii) With modest non-axisymmetric density perturba- 
tion, angular momentum transfer by gravitational torque 
can prevent fragmentation. However, even larger non- 
axisymmetry causes fragmentation through the bar-mode 
instability. 

(iv) For Z ^ 10~ 5 Zq, majority of fragmentation oc- 
cur just after the first-core formation, while for Z ^ 

Zq fragmentation epochs distributes in a wide range of 
10 14 cm~ 3 < n c < 10 22 cm~ 3 . Clouds with lower metallicity 
fragments at a higher density with shorter orbital period and 
their binary separations are narrower. Fragmentation pro- 
ceeds mainly via the bar configuration in low (< 10~ 3 Zq) 
metallicity clouds, while via the ring configuration for higher 
metallicity. The difference in fragmentation modes may af- 
fect the subsequent evolution in the accretion phase. 

(v) With smaller threshold value for the rotation param- 
eter, more binaries are expected to form from lower metal- 
licity gas. With typical rotation rate for nearby star-forming 
clouds, most stars with Z < 1Q~ 4 Zq are born as members of 
binary /multiple stellar systems. Fragments have the masses 
of 10 -3 — 0.1 Mq at their formation. If they are ejected 
from the host cloud by three-body interaction before gain- 
ing a large amount of mass, substellar-mass objects with ex- 
tremely low- or even zero-metallicity can form. The observed 
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stars with [Fe/H]<-5 might have formed by this mecha- 
nism. Moreover, we expect the possible discovery of zero- 
metallicity stars in the Milky- Way halo. On the other hand, 
if not ejected, the fragments accrete the ambient gas. In low- 
metallicity clouds, because the accret ion does not stop until 
the central stars grow to > IOOMq ijOmukai fc PallallioOll . 
120031 ; iMcKee fc Tanlliooih . massive stars or massive binary 
systems form. 
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Binary Formation with Different Metallicities 



Table 1. Fragmentation conditions for each metallicity 



Metallicity 


(n adi /10 4 cm" 3 ) 


/3 ,crit (nini/10 4 cm" 3 ) 1 / 3 




3 x 10 7 


3.2 x 10" 4 




10 s 


2.2 x icr 4 


io- 2 z Q 


7 x 10 8 


1.1 x icr 4 


io- 3 z Q 


10 9 


1CT 4 


1CT 4 Z 


10 10 


4.6 x 1CT 5 


icr 5 z Q 


10 11 


2.1 x icr 5 


1O- 6 Z 


10 13 


4.6 x 10~ 6 




n [cm 3 ] 



Figure 1. Thermal evolution of collapsing clouds with different 
metallicities (Z = 0, 10" 6 , 10 -5 , 10~ 4 , 10~ 3 , 10~ 2 , 1CT 1 and 
Zq) against the number density. To stress the variation of pres- 
sure with density, P/nk^ = T/fj. is plotted, where fi is the mean 
molecular weight. 
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10 10 1Q 12 10" 10 16 10 1B 10 [Q 10 ia 10" 10 1S 10 
r [cm] r [cm] 



Figure 2. Evolution of radially-averaged density profiles for 
clouds with different metallicity. Profiles at different times are 
plotted in the same panel for the same cloud. Formation epochs 
of the protostar and first core are indicated by solid and dotted 
arrows, respectively. The relations r~ 2,2 and r~ 2 are also shown 
in the top right and bottom left panels, respectively. 
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Z=Z 




Figure 4. Final outcomes of cloud collapse for metallicity 
Z = Zq for different combinations of the initial amplitude of non- 
axisymmctric perturbation and rotation parameter /3o- The 
density distribution (colour- scale) around the centre of the cloud 
on the equatorial plane is plotted in each panel. The grid level I 
and grid scale are shown at the upper-left and lower-right corners, 
respectively, of each panel. The colours of the panel frame indicate 
the classifications: red: fragmentation, violet: merger, blue: non- 
fragmentation, and grey: stable-core models. The white-and-black 
dotted line indicates the border between fragmentation (fragmen- 
tation and merger models) or not (non-fragmentation and stable- 
core models). 
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Figure 5. Same as Fig. [4] but for Z = 10 Zq with gray-scale. 
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Figure 6. Same as Fig.g] but for Z = 10~ 2 Z C 
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Figure 7. Same as Fig. [5] but for Z = 1CT 3 Z G 



12 M. N. Machida et al. 



10' 



m 


H 












15 * 
f 

J 

■ 1fc*J 










m r 
/ 




VI 






IliEJ 


V 


V 






• 

V 

| 

a 


"<> 

UIMJ 


r -l: a i 




2£ 

-J 

r 

kUiLiUJ 




1 






— »=w 


\ 


■ 

—-rev 


* 

— H£ 


$ 




r \ 

! V 


11 


CLZ3JJJ 




27 

■ 11 MJ 




27 

§ 

LIIAJ 


27 

■ 11 AJ 




2fi 

/ 


s 




a 

• 




ZB 

* 

tor* 


Z9 

• 


* 

*11*J 


| 

W1 *J 


§ 


?? 


1 



a.i .'.l ojl o.jl 



Ac, 



Figure 8. Same as Fig.H but for Z = 10" 4 Z G 
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Figure 9. Same as Fig. GO but for Z = 1CT 5 Z G 
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Figure 11. Same as Fig. [5] but for Z = 




Figure 12. Evolution of the normalized angular velocities at the 
centre against the number density. Models with the same (= 
0.01) but different j3 ^lO" 1 -^ 6 ) are plotted in the same pan- 
els for the same mctallicities. The cases with/3 = 10-\10- 3 and 
10 — 5 are indicated by solid lines, and those with @o = 10 — 2 , 10 -4 , 
and 10 — 6 by dashed lines. The relation oc p 1//6 , valid for spherical 
collapse, is also shown by a solid line in each panel. The verti- 
cal broken line indicates the epoch when the gas first becomes 
adiabatic. 




Figure 13. The number density and separation of fragments at 
the fragmentation epoch. All fragmentation models are plotted for 
each metallicity. The separation corresponds to the furthermost 
distance between fragments. The number of fragmentation models 
is given at the lower left in each panel, while the number of merger 
models is indicated by the figure in parentheses. The Jeans length 
is also illustrated by solid lines. The vertical broken line indicates 
the epoch when the gas first becomes adiabatic. 
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Figure 14. Histogram of the binary separation r sop at fragmen- 
tation for all fragmentation models. The metallicity is indicated 
in each panel. 



